Quadrupole-scissors modes and nonlinear mode coupling in trapped two-component 

Bose-Einstein condensates 
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We theoretically investigate quadrupolar collective excitations in two-component Bose-Einstein 
condensates and their nonlinear dynamics associated with harmonic generation and mode coupling. 
Under the Thomas-Fermi approximation and the quadratic polynomial ansatz for density fluctua- 
tions, the linear analysis of the superfluid hydrodynamic equations predicts excitation frequencies 
of three normal modes constituted from monopole and quadrupole oscillations, and those of three 
scissors modes. We obtain analytically the resonance conditions for the second harmonic generation 
in terms of the trap aspect ratio and the strength of intercomponent interaction. The numerical 
simulation of the coupled Gross-Pitaevskii equations vindicates the validity of the analytical results 
and reveals the dynamics of the second harmonic generation and nonlinear mode coupling that 
lead to nonlinear oscillations of the condensate with damping and recurrence reminiscent of the 
Fermi-Pasta-Ulam problem. 

PACS numbers: 03.75.Kk, 03.75.Mn 



I. INTRODUCTION 

The study of collective excitations gives us insights 
into the fundamental properties of many-body systems. 
The experimental observations of collective excitations 
in alkali-atomic Bose-Einstein condensates (BECs) 0, Q 
have stimulated theoretical investigations in this context 
0. When the amplitude of a collective mode is increased, 
the dynamics of a BEC enters a nonlinear regime affected 
by the interatomic interaction, giving rise to a significant 
frequency shift and mode coupling in the collective exci- 
tations. Compared with a single-component BEC, even 
richer collective dynamics have been predicted to occur in 
two-component BECs due to the intercomponent interac- 
tion [j, yj la, mM, Mi A mixture of BECs can be produced 
experimentally by simultaneously trapping atoms in dif- 
ferent hyperfine states [ToL ITll [Tj or of different species 
[lflj . Collective excitations and effects of the intercom- 
ponent interaction have been investigated experimentally 
by exciting a dipole mode, i.e., a center-of-mass oscilla- 
tion, of two-component BECs [13, HH LUj ■ The dynamics 
of the center-of-mass oscillation was studied numerically 
by Sinatra et al. they found two dynamic regimes 
of a periodic oscillation and a damped oscillation with 
strong nonlinear mixing of excited modes. The present 
paper investigates excitations of quadrupole and scissors 
modes in two-component BECs with a particular empha- 
sis on the nonlinear dynamics due to the mode coupling 
between these modes. 

The frequencies of collective modes are accurately de- 
scribed by the linear analysis of the Gross-Pitaevskii 
(GP) equation. In the Thomas-Fermi (TF) limit, which 
applies to condensates with a large number of atoms, 
the GP equation reduces to a hydrodynamic equation, 
whose analytical solution for the excitation frequencies 
can be obtained by expressing the density fluctuations 
in terms of polynomials of degree n in the Cartesian 



coordinates x 
e.g., 8n 



The polynomial of n = 1, 

,x 2 ,x 3 ) = (x,y,z), de- 



ll and z 

(1) = Y^,i=iPi x i> where (xi, 
scribes three dipole modes. The dipole modes have the 
frequencies of the harmonic trapping potential, indepen- 
dent of the interatomic interaction due to Kohn's theo- 
rem. The polynomial of n = 2, e.g., Sn 1 - 2 ^ = ^2pijXiXj, 
gives six normal modes. In an axisymmctric harmonic 
potential, linear combinations of the diagonal compo- 
nents p xx , p yy and p zz describe three normal modes: one 
m = 2 mode and two radial-breathing modes with m = 0, 
where m is the projected angular momentum on the sym- 
metry axis; we shall refer to these modes as diagonal 
quadrupole modes. The remaining three normal modes 
are associated with the off-diagonal components p xy , 
and p zx ; we shall refer to them as scissors modes |l4j . 
by which superfluidity of a BEC has been demonstrated 
|15| . The nonlinear term in the GP equation couples 
these modes more strongly for larger amplitudes of the 
excitation. The nonlinear mode coupling is also sensitive 
to the trap geometry because the excitation frequency 
depends on the trapping frequencies. Hechenblaikner et 
al. observed the second harmonic generation and mode 
coupling between the two m = radial-breathing modes 
in a single-component BEC by varying the aspect ra- 
tio of a trapping potential 0] and obtained the results 
consistent with the theoretical prediction [ij], [l8( . The 
mode couplings associated with the scissors modes have 
also been studied experimentally 0] and theoretically 

mm 

While the collective excitations in a single-component 
BEC have been well understood, much remains to be in- 
vestigated for the collective excitations of two-component 
BECs. In this paper, we first focus on the quadrupole 
and scissors modes of two-component BECs, and discuss 
their mode couplings by numerical simulations of the GP 
equations. While collective excitations of two-component 
BEC systems have been studied both analytically 
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and numerically [(J [7], |8| , no special attention has been 
paid to quadrupole and scissors modes. From a theoreti- 
cal point of view, it is not easy to predict analytically the 
collective excitation frequencies of two-component BECs 
in a trapping potential because of the complex ground- 
state structure 01111130. The assumption that 
two components overlap in the ground state allows us to 
use the polynomial ansatz for density fluctuations within 
the TF approximation and to analytically derive the fre- 
quencies of the collective modes. The quadratic polyno- 
mial ansatz gives us a clear understanding on three di- 
agonal quadrupole modes and three off-diagonal scissors 
modes. The spectra of these modes are shown to bifur- 
cate into the in-phase and out-of-phase oscillations due 
to the intercomponent interaction. Out numerical simu- 
lation of the coupled GP equation vindicates the validity 
of the analytical results. We obtain the resonance con- 
ditions of the second harmonic generation, in which two 
quanta of the fundamental mode are converted into one 
quantum of the second harmonic, depending on the trap 
aspect ratio and the strength of intercomponent inter- 
action. Our numerical simulation clarifies the dynamics 
of the second harmonic generation and nonlinear mode 
coupling. 

This paper is organized as follows. SectionlTTldescribes 
the basic equations of the system and briefly reviews 
the equilibrium properties of the system in the TF ap- 
proximation. Under the assumption of the miscibility of 
two-component BECs and the quadratic polynomials for 
the density fluctuations, the frequencies of the diagonal 
quadrupole modes and the scissors modes are found in 
Sec. IIIII In Sec. IIV Al the resonance conditions of second 
harmonic generation of the normal modes are obtained. 
The nonlinear dynamics associated with the mode cou- 
plings is numerically studied in Sec. IIVBI Section [Vj 
summarizes and concludes this paper. 



II. BASIC EQUATIONS AND EQUILIBRIUM 
STATE 



We consider two-component BECs with the atomic 
mass rrii in the hyperfine state \Fi,Mpi). The conden- 
sates are trapped by harmonic potentials which are de- 
scribed in terms of the magnetic field |B(r)| as 25] 

Vi(r) = fi B9i M Fi \B{r)\ 



pLB9iM Fi { B + -(K x x A + K y y A + K z z A 



= V? + -m^x 2 +u 2 y y 2 + u 2 z z 2 ), i = l,2, (1) 

where /is is the Bohr magneton, gi the g-factor of the i-th 
component, and LOik (k = x, y, z) the trapping frequency 
satisfying the relation 



m l J 2 k = fj, B giM Fi K k 



(2) 



In this paper, we do not consider the relative displace- 
ment of the potential minima caused by the constant 



term V® — (XBgiM Fi Bo by assuming g\Mp\ = g 2 M F2 
[26|. We also ne glec t the displacement of the system 
due to gravitation [jjj] • The dynamics of two-component 
BECs are described by the coupled GP equations for the 
condensate wave functions ipi and 2 , 





' h 2 v 2 


lh -dT = \ 


v, 2toi 




f n 2 v 2 


lh -df = { 


V 2m 2 



+ Vi+ M i|0i| 2 +ui 2 |V> 2 | 2 JV'i, (3a) 
+ V 2 + u 2 \iJ2\ 2 + u^lV'il 2 ) 4>2. (3b) 



Here the intracomponent interactions (iti and u 2 ) and 
the intercomponent one (ui 2 ) are expressed in terms of 
the corresponding s-wave scattering lengths oi, a 2 and 
ai 2 as Ui = Anh 2 di / mi (i = 1,2) and u\ 2 — 2nh 2 a\ 2 j 'm\ 2 
with m\ 2 = mim 2 /(mi + m 2 ) being the reduced mass. 
In this work, we will discuss the case of positive scat- 
tering lengths. By writing the condensate wave function 
as ipi(r, t) = yjrii{r, i)e* 9i ( r '*), the time-dependent equa- 
tions for the condensate density n.j(r, t) and the velocity 
field Vj = {%/mi)V9i are given from Eqs. (J^i) and (J2t>) 
as 



dni 



-v-Cmvi), 



mi- 



9vi 
~dt 



,_. h 2 v 2 Vm Tr 

-VI Z=- + Vi-m 



2mi y / nT 



1 



+u\ri\ + u 12 n 2 + -miv 1 ) , 



m 2 - 



dv 2 
~dt 



= -V 



dn 2 
~dt 
h 2 V 2 y/m 
2m 2 Jn^ 



V ■ (n 2 v 2 ), 



V 2 - fj,2 



+u 2 n 2 + u 12 ni + -m 2 v 2 



(4a) 

(4b) 
(4c) 

(4d) 



The chemical potential /ij is determined from the nor- 
malization condition Ni = J drrii where Ni is the number 
of particle of the i-th component. 

The frequencies of collective modes are calculated by 
linearizing Eqs. around the equilibrium solutions as 

rii = ni + 5ni and v, = 8vi . The equilibrium solutions 
of the cou pled GP equations show a variety of stable 
structures 1IJ El El El EI When the intercompo- 
nent interaction is strongly repulsive, the two compo- 
nents separate spatially; the condition of the phase sepa- 
ration for trapped two-component BECs can be derived 
simply with the TF approximation |2ll |24, 25] , in which 
the quantum pressure terms (?i 2 /2miy / rI7)V 2 y / n7 in Eqs. 
|J1]d) and J2JI) are neglected. In the spatial region where 

two components overlap (nf^ ^ and n!, ' ^ 0), the 
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equilibrium density distributions are given by 



effective when the TF radii 



(0) 



U1U2 — u 



12 



M2-r 2 Vi(r) 

Miri2 L u 2 



UlU2 — U 



1.2 



u 2 r 



1.2 



U\2 



l H-T 1 V 2 {r) 



(5a) 



(5b) 



Here we dehned the following dimensionless parameters, 



r X2 = 1 - 



mi + m 2 2 

4mi2 aia 2 ' 
mi a 12 



1 



r 2 = i- 



2m i2 ai 
w 2 ai2 

2mi2 a 2 



(6) 
(7) 
(8) 



In order for the two components to overlap, the right- 
hand sides of Eqs. J5]i) and (J5j}) must be positive. The 
overlapping region of the two components is delimited 
by the boundary defined as fa — u\2^j/uj — TjVi(r) = 



In addition, the parameter Ti 2 must be positive, 
for if Ti2 < 0, i.e., U1U2 — u\ 2 < 0, there is no overlap- 
ping region, which implies the phase separation of the 
two components. The phase -sep arated solutions will not 
be considered in this work [23. For Ti2 > 0, there are 
two possible overlapping density profiles depending on 
the signs of the parameters Ti and T 2 . These signs de- 
termine the curvature of the density profile at the center 
of the trapping potential because Ti and T 2 are the co- 
efficients of the x 2 , y 2 and z 2 terms in Eqs. (jHfc.) an d 
JHb)- For Ti > and T 2 > 0, both density profiles are 
written by inverted parabolas in the overlapping region. 
For Ti > and r 2 < (T x < and T 2 > 0), one density 
profile is an inverted parabola, but the other is a usual 
parabola in the overlapping region. In the peripheral re- 
gion in which the density of one component vanishes, the 

other component follows the single-component TF pro- 
Jo) 



pTF _ 



' Mi - U12H2/U1 
T 2 miUj 2 k /2 



pTF 
H 2k 



\ [12 - U 12 fll/u 2 

r 1 m 2 w| i ./2 



(9) 

(fc = x, y, z) of the two overlapping condensates are al- 
most equal: 



pTF 
H 2k 
RTF 



U2 r 2 n 2 



1/5 



1. 



(10) 



In addition, the amplitude of the fluctuation should be 
small for the assumption to be valid. As shown later, 
the predicted frequencies of collective modes are indepen- 
dent of the particle number under the TF approximation. 
Therefore, for any and ai the condition (|ll)fl can be 
met by adjusting the ratio of the particle number N2/N1. 

The frequencies of collective modes are obtained by 
linearizing Eqs. Qi) and (QJs) around the equilibrium 
solutions as rii = nf^ +6rii and Vj — <5v^. Here we neglect 
the quantum pressure term in Eqs l@J. This results from 
the assumption that the spatial variations of the density 
are smooth not only in the ground state but also in the 
low-lying excited states. The linearized equations then 
read 



d 2 S ni 



1 



- — V-{r4 0) V(ui<5ni+wi2<5n 2 )}, (Ha) 
ot mi 

^Ji = _L V ■ {n^V(u 2 5n 2 + u 12 Sm)}. (lib) 

Ot z TO 2 

The analytic solutions to these equations are obtained by 
expressing the density fluctuations in terms of quadratic 
polynomials 

Snx(r,t) = p {t) -p x (t)x 2 -p v (t)y 2 - p z (t)z 2 

-Pxy(t)xy - p yz (t)yz - p zx (t)zx, (12a) 
Sn 2 {r, t) = q (t) - q x {t)x 2 - q v {t)y 2 - q z {t)z 2 

-q xy (t)xy ~ q yz (t)yz - q zx (t)zx. (12b) 



file m = {fii - Vi)/Ui. This profile is connected to the The diagonal qua drupole modes are described by lin- 
over lapping density profile of Eq. JoJ by adjusting the 
chemical potential (j+ (i = 1 or 2) under the normaliza- 
tion condition J dvnf = Ni [2l| . 



ear combinations of the diagonal components p x ,y,z and 



III. FREQUENCIES OF THE LOW-LYING 
COLLECTIVE MODES 



In this section we discuss the excitation frequencies of 
the diagonal quadrupole modes and the scissors modes 
of the two-component BECs by solving Eqs. Q}. In the 
following analysis, we assume that the two components 
overlap completely, use Eqs. (0i) and (03) with inverted 
parabolas (ri > and T 2 > 0) as the equilibrium solu- 
tion , and apply no boundary conditions to the so- 
lutions of Eqs. 10^) and (0J}). This assumption may be 



y,z 

1x,y,z, while the off-diagonal components ,,,, and 
qxy,yz,zx determine the scissors modes [lil Il8l l2Qj . The 
time dependences of po and go are determined by the con- 
ditions J dvSrii = 0, which express the particle-number 
conservation of each component. The equilibrium states 
n[ ^ are given by Eqs. JSJl) (JSJs). It follows then from Eqs. 
(|llb ) and (|llb ) that the time-dependent variables X = 

{Px,q x , Py , qy , Pz, qzi Pxy t qxy, Pyz j qyziPz x ,q zx ) T obey 



dt 2 



X = -M-X = - 



M Q 
M s 



•X. 



(13) 



In the linear limit, the diagonal quadrupole modes and 
the off-diagonal scissors modes are completely decoupled 
IHI23. The matrix M Q is too large to be written down 
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here, but Ms for the scissors modes is simply given as 



M xy 





! 













(14) 










M zx J 







with 

M 3k = 



Ti U12 I, ,2 , , ,2 \ 




.(15) 



Since the three scissors modes are also decoupled as seen 
from Eq. (|14f> . we can treat the xy scissors mode sepa- 
rately. The frequency Sl xy is easily calculated from M. xy 



by substituting p xy = Tie 1 



and q x 



T 2 e i: 



into 



Eq. (| 1 31) and solving the resulting secular equation with 
the result 



n 2 =-L 
xy 2r 15 



r 2 (a4 + <) + r 1 (u4 + <4J 



±4 



[r 2 (w? a 



^ ri r 2 (^ + O^L + ^) 

1/2-1 



■ w 



Ti(w 



(16) 



For the case of equal atomic mass m\ = m 2 (hence oj\k 
u>2k = by Eq. (J2J), Eq. I|16[l is simplified to 



fl xy+ = Jul + u. 



^xy— 



■ w 2 v 5a 



with 



5a 



{a% - ai 2 ){a 2 - a 12 ) 



aia 2 



(17a) 
(17b) 

(18) 



When the two components overlap completely, we have 
< 5a < 1, so that £l xy - < Vl xy+ . Thus, the frequencies 
of the scissors modes split into two branches. The up- 
per frequency £l X y+ correspond to the in-phase oscillation 
which keeps the principal axes of the two components in 
line. For mi — m 2 , the frequency of the in-phase mode is 
independent of the intercomponent interactions because 
of Kohn's theorem 0. The lower frequency Q X y- corre- 
sponds to the out-of-phase motion of the two components 
in which they oscillate with the 180° phase difference. 
The frequency of the out-of-phase mode depends strongly 
on the interatomic interaction; the measurement of both 
frequencies of the in-phase and out-of-phase modes thus 
allows us to determine 5a from Eqs. (|17b ) and (117b ). 

The frequencies of the diagonal quadrupole modes are 
also calculated by solving the corresponding secular equa- 
tion of Mq. For m\ = m 2 and the axisymmetry of a 



trapping potential lu x = uj y / X 
for Mq yields 



n 2 + 
n 2 _ 

OJ 2 



the secular equation 



(19a) 



16A 2 + 16 ' 



n 2 



-5a. 



(19b) 




FIG. 1: Schematic representation of the quadrupole and scis- 
sors modes of two-component BECs. Oscillations of compo- 
nent 1 are indicated by solid arrows and those of component 
2 by dashed arrows. Mode coupling processes in Table Q are 
also shown. 



First, we quote the excitation of a single component 
condensate, i.e., ai 2 = (5a = 1). When A = \/8 and 
Qi2 = 0, for example, Eqs. fl^|l give Q/uj x = \/2, 1-79 
and 4.97 [H|, which corresponds to a quadrupole-type 
oscillation in the x-z plane with m = 2, an in-phase os- 
cillation in the x and z directions but out-of-phase in 
the y direction, i.e., the m = low-lying mode, and an 
in-phase compressional mode in all directions i.e., the 
m = high-lying or breathing mode, respectively. Three 
scissors-mode frequencies are given by il zx /u> x = \/2 and 
^xy/^x = Qyz/u x — 3; the zz-scissors mode is degener- 
ate with the m = 2 quadrupole mode . 

For ai2 the frequencies of all modes split into two 
branches as mentioned above; the corresponding modes 
are graphically represented in Fig. In Fig. |2fa), we 
show the frequencies of twelve collective modes as func- 
tions of 012/01 for the condensates with a\ = a 2 . The 
frequencies of the out-of-phase modes decrease monoton- 
ically with ai 2 , because the stronger mutual repulsion fa- 
vors the phase separation and makes these modes softer. 
The frequencies vanish at 012/01 = 1, at which the lin- 
ear analysis around n'p of Eqs. and JSJd) with an 
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FIG. 2: ais/ai dependences of the twelve collective normal 
modes obtained from Eq. 1131 for the parameter ai = d2, 
A = \/8, and (a) m,2 = mi and (b) — 0.75mi. Each mode 
splits into the in-phase mode (thin curves) and the out-of- 
phase mode (dashed curves) for ai2 7^ 0. In (a), the numeri- 
cally obtained frequencies are plotted by open symbols (for in- 
phase) and filled symbols (for out-of -phase) . As explained in 
Sec. IIV Bl the reduction to a two-dimensional system changes 
the frequencies of the m = low- and high-lying modes; for 
A = Vo their values become small from those of Eqs. 1191 by 
10 % for the low-lying mode and 1% for the high-lying mode. 

inverted parabola breaks down. 

For mi 7^ m 2 the excitation frequencies of all modes 
separate even at a\i = 0, because the trapping frequency 
of one component is different from that of the other as 
seen in Eq. @. The frequencies for m 2 = 0.75m! are 
shown in Fig. 0b). In this case, the frequencies corre- 
sponding to the in-phase motion decrease gradually with 
a\2, but the qualitative behavior is nearly the same as 
that in Fig. Eta). In general, the frequency of the out- 
of-phase modes vanishes when 

rir 2 = 1--^— 1--^— =0, 20 

\ 2mi2 at J \ 2m 12 a 2 J 

which gives, e.g., a^/ai = 0.856 for m 2 = 0.75mi. These 
results are consistent with the numerical study of the 
excitation frequency Q , where the frequency of the out- 
of-phase mode also decreases with ai 2 and the ground 



state undergoes phase separation eventually. 

IV. HARMONIC GENERATION AND MODE 
COUPLNG AMONG QUADRUPOLAR MODES 

In the preceding section, we obtained the frequencies 
of the diagonal quadrupole modes and the scissors modes 
with the in-phase or out-of-phase motions of two compo- 
nents. Although these modes are decoupled in the linear 
limit, the growth of the excitation amplitude makes the 
nonlinear coupling among these modes stronger, casting 
the dynamics into a nonlinear regime. The mode cou- 
pling also depends strongly on the geometry of a trapping 
potential (la. XVA Ha . il9l . |20| . In this section, we will show 
the nonlinear dynamics caused by the harmonic genera- 
tion and the mode coupling among the quadrupole exci- 
tations by numerically solving the coupled GP equations 
The nonlinear term in the GP equation is propor- 
tional to |-0| 2 , being able to generate an excitation with 
a frequency of the integral multiple of that of an ini- 
tially excited mode. This fact allows the spontaneous 
up-conversion process, while the down-conversion can- 
not be explained within the framework of GP equation 
[Til flU, 120] . In this work, we discuss the up-conversion 
process of the harmonic generation and confine ourselves 
to the second-harmonic process, where the resonance 
of two modes occurs via annihilation (creation) of two 
quanta in a lower-energy mode and creation (annihila- 
tion) of one quantum in a higher-energy mode i.e., sec- 
ond harmonic generation at l^hi p-h = 2£li nw . Compared 
with the single-component case |l6l Int Il8l Il9l |2C| , an 
even richer variety of mode couplings can occur because 
of the presence of the out-of-phase modes, caused by the 
intcrcomponent interaction. 

A. Resonance condition 

Before discussing the numerical results, we evaluate 
the resonance condition between two quadrupole modes 
derived in the preceding section. In order to simplify the 
problem, we consider two-component condensates with 
the equal mass (mi = m 2 and oj\k = W2fc), trapped in an 
axially symmetric potential (lo x = <~o y j\ = lo z ). Then, 
simple algebraic manipulations of Eqs. i|17|) an d ESt 
yield some resonance conditions for the mode coupling. 

1. Coupling between two diagonal quadrupole modes 

First, we consider the mode coupling between two di- 
agonal quadupole modes. To take account of the non- 
linear effect, we follow the method of Ref. 0]; the 
density and the velocity are assumed to take the forms 
ni(r,i) = po{t)-p x (t)x 2 -p y (t)y 2 -p z (t)z 2 andvi(r,t) = 
V{a x (t)x 2 + a y (t)y 2 + a z (t)z 2 }, and similar forms for 
component 2 by replacing p — > q and a — > /3, where we 
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neglect the off-diagonal contribution. The substitution 
of these terms into Eqs. (QJ yields the equations of mo- 
tion for the time-dependent variables p, q, a and /?. As 
in Sec. IIIII we impose no boundary condition on these 
variables. These equations are simplified by introducing 
the new variables ak and bk (k = x, y, z) defined by ^3] 



(a) 10 



Pk 



2u 1 Ti 2 a x a y a z a\ ' * 2u 2 Ti 2 b x b y b z b\' 
Then we obtain otk = hk/2ak, (3k — bk/2bk, and 



Y x m 2 ijj\ k 



a-k 



-2k 



i 



-Clk 



-b k 



"121"! 



T 12 a x a y a z a k u 2 T 12 b x b y b z b 2 k 
Ti 1 , u i^2 b k 
r 12 b x b y b z b k uiTi 2 a x a y a z al 



(21) 

,(22) 
(23) 



The linearization of Eqs. IL'l'l) and (|23|) yields the matrix 
Mq, whose eigenvalues are the same as those of Mq 
of Eq. H13J) . In order to clarify the mode coupling, we 
rewrite ak and bk by using the new basis obtained by 
the eigenvectors e, of Mq as x = $^ i=1 Cj(i)e,*, where 
x = {a x ,a y ,a z ,b x ,b y ,b z ). From Eqs. i(2"2)l and i(2"3|) one 
obtains the equations for the time-dependent variables 
Ci(t); these equations are completely decoupled in the 
linear limit. 

In order to understand the process of second harmonic 
generation, we expand the equations of motion up to sec- 
ond order in Ci(t), obtaining the form 



-fife* 



^IjkCjCk 
3,k 



(24) 



with the mode frequency f2j. The second- harmonic pro- 
cess between the i- and j-modes occurs via the terms like 
JjjCj (i j) for Qi ~ 2flj. For a single-component con- 
densate (ai2 = 0) in an axially symmetric potential, one 
can easily show that at A = yl6/7 the second harmonic 
of the to = 2 mode coincides with the m = high-lying 
mode, and at A = (5v / 5±V / 29)/6v / 2 the second harmonic 
of the to — low-lying mode coincides with the m = 
high- lying mode |17| . For two-component condensates 
with mi = TYi 2 , these conditions for A also hold for the 
mode couplings between in-phase modes. 

We summarize the other possible second-harmonic pro- 
cesses in Table H] and Fig. ^ The condition of the mode 
coupling is obtained by Q$ = 2Qj in Eqs. l(T§|l . being 
represented by the aspect ratio A and the interaction pa- 
rameter Sa as shown in the middle column of Tableland 
in Fig. 01 a). Since A > and < Sa < 1, depending 
on the relation between two frequencies, each coupling 
process occurs within a certain range of Sa as shown in 
the right column in Table [I] The processes (III) and 
(VI), which describes the coupling between the in-phase 
and out-of-phase modes with m = mode, occur for ar- 
bitrary A when Sa — 0.25. For condensates with equal 
mass there is no up-conversion process in which higher- 
lying modes are the m — 2 modes and the out-of-phase 
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FIG. 3: Conditions of the second-harmonic resonance in the 
parameter space of A and Sa = (ai — avi){a.2 — ai2)/(aiia22 — 
a\ 2 ). The mode couplings between two diagonal quadrupole 
modes are shown in (a), and those between a diagonal 
quadrupole mode and a scissors mode in (b). 



to = modes, because the term jjjCj does not appear 
for these modes. (For different mass mi ^ m 2l the out- 
of-phase to = modes can become higher-lying modes of 
the up-conversion process.) 



Coupling between a diagonal quadrupole 
scissors mode 



mode and a 



Now we discuss the mode coupling related with the off- 
diagonal scissors modes. When the off-diagonal contribu- 
tion is included in the fluctuation, the analysis in the pre- 
ceding subsection cannot be used because an appropriate 
transformation like Eq. (|21|l does not exist. However, if 
we expand the condensate wave function as i^i(r, t) = 
e -*w*/&[^ 9 ( r )+ lti ( r ) e -*n<* + v * (r)e iQit ], the presence of 
the mode coupling is understood by the symmetry of the 
mode functions Uj(r) and Wj(r). The coupling coefficient 
for the second-harmonic process is written by the overlap 
integral of the ground state wave function and the mode 
functions as / c?r?/; i g(r)u 1 * ow (r) 2 Uhigh(r) HEH^I- The scis- 
sors modes are characterized by the odd-parity function 
Mj(r) tx xy,yz,zx. Thus, there is no coupling between 
a higher-lying scissors mode and a lower-lying diagonal 
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TABLE I: List of all possible second-harmonic processes for two-component BECs with mi = mi in axisymmetry potentials. 
IP and OP stands for "in-phase" and "out-of-phase" , respectively. 



second-harmonic process 



aspect ratio A = uj y /ui x 



range of Sa 



(I) m = 2 or jza;-scissors, OP — > m = - low, IP 

(II) m — 2 or jz:E-scissors, OP — > m = - high, IP 

(III) m = 0-high, OP -> m = - high, IP 

(IV) m = - high, OP -> m = - low, IP 

(V) m = - low, OP -> m = - high, IP 

(VI) m = 0-low, OP -* m = - low, IP 

(VII) xy, j/z-scissors, OP — * m — 2, IP 

(VIII) xy, yjz-scissors, OP — *■ m — - low, IP 
(IX) xy, 2/2-scissors, OP — ► m — - high, IP 



<5a(l-2<5a) 
5— 125a 



<5a(l-2<5a) 



5-12<5a 

A = arbitrary 

\/ 805a 2 -8<Sa + 5±v^(l+4<5a)^/5-56<Sa+80«a 2 

6\/ <5a 

\J 80(5a 2 -8ia + 5±V^(l + 4«a)^/5-56<5a+80(5a 2 
6 v^5<T 

A = arbitrary 
A = 



16<5a 2 -14<5a+5±\/l64<5a 2 -140«a+25 
5a(3-45a) 



< Sa < 5/12 
1/2 < Sa < 1 
<5a = 0.25 
< 5a < I=2& 

i±m <5a< i 



Sa = 0.25 

< Sa < | 



< 5a < 



5(7-2^2) 



, _ 1 / 163n 2 -14an+5±Vl64<5n 2 -140i5n+25 5(7+2^2) c . J § for(+) 
2 y Sa(3-4,5a) 82 < 0a < | j for /_) 



quadrupole mode which has the even parity, because the 
overlap integral vanishes. However, there can be second- 
harmonic processes from a lower-lying scissors mode to a 
higher-lying diagonal quadrupole mode. 

The resonance condition for the zs-scissors mode is the 
same as that of the diagonal m = 2 quadrupole mode, 
because they are degenerate in an axisymmetric poten- 
tial. Thus, we focus on the mode coupling between the 
xy- and yz-scissors modes with the same frequency. For 
their in-phase mode, no resonance occurs for the diago- 
nal in-phase or out-of-phase quadrupole modes, because 
the condition fi^ = 2fij cannot be satisfied in the case 
of an axisymmetric trap. On the other hand, there are 
couplings for their out-of-phase modes as shown in (VII) , 
(VIII), and (IX) in Table [fl Fig. [J and Fig. CJb). 



3. Coupling between scissors modes 

A similar symmetry argument shows that there is no 
direct coupling via a second-harmonic process between 
two scissors modes, because the overlap integral of the 
second harmonic generation is always zero. Actually, if 
we neglect the diagonal modes, substitutions of the form 
5ni = p (t) - p xy (t)xy - p yz (t)yz - p zx {t)zx and 5a\ = 
V{— a xy (t)xy — a yz (t)yz — a zx (t)zx} etc. into Eqs. (@J 
yields 



Pxy — 


-2(pi 0) - 


^Py 0) )axy 


PzxQ-yz 


Pyz^-zx i 


(25a) 


Pyz = 


-2(H 0) - 




Pxy^zx 


Pzx^xy ■> 


(25b) 


Pzx = 




>rP { x ) )a zx 


Pyz^xy 


PxyQyz ■> 


(25c) 




&xy 


Ui 

Pxy 


Ml2 
T Qxy 




(25d) 




mi 


mi 






&yz 


Ml 

= Pyz 


H Qyz 


^zx^xyi 


(25e) 




mi 


mi 






a zx 


Ml 

= Pzx 


H q» x 


O^xy^-yzi 


(25f) 






mi 


mi 





and the similar forms by replacing all indices 1 by 2 and 
(p, a) by (q, f3), where p x °} y , z and q x °} y , z are the coefficients 
of the equilibrium condensate obtained by Eqs. (|5ji) and 
(JSJd) . These equations show that all scissors modes should 
initially be excited to induce the coupling between scis- 
sors modes. This conclusion was also obtained by Al 
Kawaja and Stoof for a single-component [2(i|. 



B. Numerical simulation 

To complement the analytical results obtained above, 
we study the dynamics of two-component BECs by the 
numerical simulations of the time-dependent GP equa- 
tions (|2J). For simplicity, we study the case of ui = a 2 
and mi = ttt.2, and let an and A = uj y /uj x be the variable 
parameters. We confine ourselves to the two-dimensional 
system, which is adequate for studying the mode cou- 
pling. A two-dimensional TF analysis of two-component 
BECs yields six normal modes; when an = 0, these are 
a a;y-scissors mode with the frequency £l xy /uj x — 1 + A 2 
and two diagonal quadrupole modes with the frequency 
Sl%/u>l = (3 + 3 A 2 T V9A 4 -14A 2 + 9)/2, correspond- 
ing to the in — low- and high-lying modes in the 
three-dimensional case. The frequencies of the out-of- 
phase modes are multiplied by a factor Sa given in Eq. 
(jT%)l . decreasing with oi 2 . The numerical calculations 
are done using an alternating direction implicit method 
as explained in Ref. |28j . Initially, we search the equi- 
librium solution of two-component BECs by the norm- 
conserving imaginary time propagation of Eqs. (JHt) & n d 
(J3t>) - Starting from this initial state with appropriate 
modulation of the phase profile, we calculate the time 
evolution of the condensates. The parameter values are 
followed by the condition in Ref. the number of 

the particles is 10 4 and the aspect ratio of the trapping 
frequency is u> y /oj x — y/8. 
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FIG. 4: Time development of (a) (x 2 )i,2 and (b) (y 2 )i,2, show- 
ing the mode-coupling process of (VI). The parameters used 
are aiz/ai = 0.665, A = yo, and modulation factor e = 0.05. 



We first numerically check the excitation frequencies 
obtained analytically in the TF approximation. In or- 
der to excite the collective modes, we change the phase 
0i{x,y) of the obtained initial state by a small factor 
ej = 0.005; for the rcy-scissors mode 9i(x,y) = eixy, 
and for the m = low- and high-lying mode 0i{x,y) — 
ylx 2 + q ± y 2 ) with q± = (-3 + 3A 2 ± V9 - 14A 2 + 9A 4 )/2 
[l7| . For the in- phase (out-of-phase) motion we set 
£i = (— )e2- We calculate the time evolution of the mean 
values (x 2 )i i 2, (y 2 )i,2, and (xy)i,2, fitting a sinusoid to 
the calculated data in order to extract frequencies, where 
the first two values characterize the diagonal quadrupole 
modes and the last one does the scissors mode. By vary- 
ing intercomponcnt interaction a\2, we plot the frequen- 
cies of the three collective oscillations in Fig. The 
numerical results agree fairly well with the analytical re- 
sults, but as ai2/ai approaches unity, the frequencies of 
the out-of-phase motion shift upward from the analytical 
ones. This is associated with the large amplitude oscil- 
lation induced spontaneously because of the softening of 
the frequencies or the deviation from the approximation 
with the overlapping TF profile. 

Next, we investigate the dynamics of the mode cou- 
pling when the second harmonic of a lower mode coin- 
cides with the frequency of a higher mode. Numerical 
simulation shows that the oscillation indeed exhibits sig- 
nificant nonlinear behavior because of the second har- 
monic generation. We show two examples of this nonlin- 
ear dynamics in Figs. 01and[SJ Figure 0] shows the time 
development of the mean widths of the condensates for 



A = \/8 and Sa = 0.20 (axi/ai — 0.665), when we excite 
initially the m = low-lying out-of-phase mode. The 
condensate widths (a; 2 )i,2 and (y 2 )i,2 undergo an out- 
of-phase oscillation, and (a; 2 )i ((y 2 )i) and (a; 2 )2 {{y 2 )2) 
undergo an out-of-phase oscillation, too. As time passes, 
keeping (x 2 )\^ and (y 2 )i.2 out-of-phase, the oscillations 
of (x 2 )i and (x 2 )2 as well as {y 2 )i and (y 2 )2 change from 
out-of-phase to in-phase as shown in Fig. 0] This time 
development shows directly the process of (VI) in Ta- 
ble [I] When we calculate under the condition Sa = 0.25 
(ai2/ai = 0.60) shown in Tabled (VI), this nonlinear be- 
havior is weaker than that in Fig. 21 The deviation of 
the resonance from Sa = 0.25 to Sa = 0.2 is due to the 
frequency shift appearing in the present numerical simu- 
lation. After a long time, the oscillation gets back to the 
original out-of-phase behavior. This is nothing but the 
recurrence phenomena known as the Fermi-Pasta-Ulam 
problem [29j; the similar behavior is also found in the 
dynamics of a rotating BEC |28|. 

The mode coupling from the out-of-phase xy-scissors 
mode to the m — high- lying in-phase mode ((IX) in 
Table UJ) is shown in Fig. [5] Figure Ufa) shows the 
time development of the mean values for component 1 
for A = and Sa = 0.613 (oi 2 /ffli = 0.24). From an 
enlargement view shown in Fig. EJb), we see that the os- 
cillation of the mean value (xy) i is damped quickly, con- 
verting into the oscillations of the diagonal component 
(x 2 )i and (y )i, where they undergo an in-phase oscil- 
lation; those of the component 2 oscillate in the same 
manner. Again, the oscillation recovers the out-of-phase 
of the xy-scissors mode. 



V. SUMMARY AND CONCLUSION 

In summary, we have studied the diagonal quadrupole 
modes and the scissors modes of two-component BECs, 
and the nonlinear dynamics associated with the sec- 
ond harmonic generation of the quadrupole and scissors 
modes. Within the TF approximation and quadratic 
polynomial ansatz for the density fluctuation, the linear 
analysis leads to analytical expressions for the collective- 
mode frequencies. These analytical results are confirmed 
by the numerical simulation of the GP equations. In- 
cluding the nonlinear contribution of the fluctuation, we 
have clarified a variety of the mode couplings between 
collective modes. The condition of the second harmonic 
generation is obtained as functions of the aspect ratio 
of a trapping potential and intercomponent interaction 
strength. At the resonance of two collective modes, the 
collective dynamics of two-component condensates shows 
remarkable nonlinear behavior involving damping and re- 
currence of the density oscillation. 

In this paper, we have confined ourselves to the con- 
densates with their density profiles completely overlap- 
ping. Two-component BECs consisting of the same 
atoms realized so far are mixtures of the states |1, — 1) 
and |2, 1) of 87 Rb in JILA [nj, |2, 1) and |2,2) of 87 Rb 
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FIG. 5: (a) Time development of {xy)i, (x 2 )i — (x 2 )^ and 
(y 2 ) i — (y 2 ) j ' , where () represents the initial value, showing 
the mode-coupling process of (IX). The parameter values are 
«i2/«n = 0.240, A = y/E, and modulation factor e = 0.2. (b) 
An enlarged view of (a) for the time interval [0,40]. 



in LENS [13 and 1 1 , 0) and |1,1) of 23 Na in MIT [T^. 
Although the scattering lengths of the mixture in JILA 
satisfy the overlapping condition, the resulting double 
condensates do not form the inverted parabolic density 
profile. The latter two mixtures lie slightly in the phase- 
separated region. Moreover, a mixture of 41 K and 87 Rb 
BECs also lies deeply in a phase-separate region . By 
changing the scattering lengths via the Feshbach reso- 
nance, however, we can, in principle, prepare the con- 
densates in an overlapping regime and verify our predic- 
tion We believe that this work will help the under- 
standing of the complex dynamical behavior of mutually 
interacting multicomponent BECs. 
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